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f^* ) We propose a model based on a generalized effective Hamiltonian for studying the effect of 

, noise in quantum computations. The system-environment interactions are taken into account 

O ' by including stochastic fluctuating terms in the system Hamiltonian. Treating these fluctu- 

ations as Gaussian Markov processes with zero mean and delta function correlation times, 
\ we derive an exact equation of motion describing the dissipative dynamics for a system of 

I , n qubits. We then apply this model to study the effect of noise on the quantum telepor- 

■ tation and a generic quantum controlled-NOT (CNOT) gate. For quantum teleportation, 

the effect of noise in the quantum channels are found to be additive, and the teleportation 
C^l \ fidelity depends on the state of the teleported qubit. The effect of collective decoherence 

i — i ! 

, is also studied for the two-qubit entangled states. For the quantum CNOT gate, we study 

the effect of noise on a set of one- and two-qubit quantum gates, and show that the results 
can be assembled together to investigate the quality of a quantum CNOT gate operation, 
■y | We compute the averaged gate fidelity and gate purity for the quantum CNOT gate, and 

investigate phase, bit-flip, and flip-flop errors during the CNOT gate operation. The effects 
of direct inter-qubit coupling and fluctuations on the control fields are also studied. We find 



that the quality of the CNOT gate operation is sensitive to the strengths of the control fields 
and the strengths of the noise, and the effect of noise is additive regardless of its origin. 
We discuss the limitations and possible extensions of this model. In sum, we demonstrate a 
simple model that enables us to investigate the effect of noise in arbitrary quantum circuits 
under realistic device conditions. 



I. INTRODUCTION 

Quantum information processing is of much current interest^]. The realization of quantum 
algorithms using nuclear magnetic resonance (NMR) DDD3 and ion-trap (J techniques has 
shown that quantum computing is possible in principle. Recent efforts for building quantum com- 
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puters have focused on techniques based on solid-state devices that are believed to be more scalable 

m 

. However, such solid-state devices usually require sophisticated manufacturing techniques, 
and the inevitable interactions between a qubit and its surrounding environment ("bath") intro- 
duce noise into the quantum system, resulting in the degradation of the quantum superposition 
state. Thus, the extra degrees of freedom of a solid-state system and the inherent system-bath 
interactions pose a great problem for quantum computing with such devices. The decoherence 
problem is the main obstacle towards the realization of a universal quantum computer, and a 
sound theoretical framework for the description of the decoherence and population relaxation of 
qubit systems is necessary Q] . 

Because the ability to compute and predict the behavior of a quantum circuit under the influence 
of noise is crucial, a model that can describe errors from the system-bath interactions could be 
extremely useful. Such a model will also be useful in the study of quantum error-correcting and 
error-preventing schemes, as well as provide informative guidelines for the design of quantum 
computers. However, describing the non-equilibrium decoherence and population relaxation of a 
many-qubit system is non-trivial. No general model exists for this purpose. Classical noise models 
and microscopic noise models have yielded some success, but these formulations do not provide a 
general solution framework for a many-qubit system. 

Classical noise models that describe the decoherence and population relaxation as exponential 
decays of the off-diagonal and diagonal components of the density matrix are widely used for the 
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3, 



but generally these models lack 



estimate of the error rates during quantum computation 
quantum features that are important for quantum computing, such as the quantum interference 
effect. 

Microscopic noise models based on the spin-boson Hamiltonian that explicitly include the linear 
couplings between t he sy stem and the bath degrees of freedom have provided valuable insights about 
decoherence effects |lfll. Illl. llfll). Recently, the decoherence and gate performance of a quantum 
controlled-NOT (CNOT) gate operation for several different physical realizations has been studied 



based on such spin-boson type Hamiltonians 

been developed to solve dynamics of microscopic Hamiltonians 3|. However, these methods are 
often complicated, and difficult to generalize for systems with more than two qubits. In addition, 
in many cases the exact form of the system-bath interactions is unknown, or the parameters are 
difficult to obtain experimentally, and the microscopic models are difficult to use in these cases. 

The Bloch-Redfield formalism is generally used to study NMR spin-dynamics ^{J, and has been 
applied to study the dynamics of many-spin systems 20]. However, this formalism, while suitable 
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in NMR systems, is not always applicable in general qubit systems. Moreover, the Bloch-Redfield 
formalism is also known to violate the complete positivity of the reduced density operator at short 
times. To apply the Bloch-Redfield formalism to quantum computing, non-physical additional time 
intervals have to be inserted between the switching events 1J, |15|]. These extra time periods will 
result in the over-estimation of the degradation of the qubit systems. 

Thus, generally speaking, a method that can be used to analyze the quality of a functional 
quantum circuit and capable of providing a quantitative result is still not available. In this work, 
we propose a stochastic Liouville equation approach to describe errors in quantum computations. 
This approach originates from the Haken-Strobl-Reineker (HSR) model first proposed by Haken 
and Strobl and later extended by Reineker in the 1970s to describe charge and energy transfer 
in organic crystals 
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|22j, 23]. own to be able to capture the coherent and 

incoherent dynamics of quantum two-level systems. In this model, the system-bath interactions are 
taken into account by allowing the site energies and the off-diagonal matrix elements of the system 
to fluctuate over time. We generalize the idea of Haken and Strobl to describe a system of n qubits. 
The resulting stochastic Liouville equation is then solved to obtain a set of equations describing the 
dynamics of a general n qubit system. To demonstrate the applicability of our method, we study 
the effect of noise on quantum teleportation and a generic CNOT gate operation, and then compare 
our results with previous work. We show that our model can reproduce the main results obtained 
previously by using microscopic model Hamiltonians. The limitations and possible extensions of 
our semiclassical model are also discussed. 



II. THE STOCHASTIC LIOUVILLE EEQUATION APPROACH 



Previous work on the study of the population relaxation and decoherence of qubit systems is 
usually based on the spin-boson Hamiltonian, in which the qubits are coupled linearly to the bath 



degrees of freedom (t 



oscillators 
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le environment), and the bath is treated explicitly as a system of harmonic 
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]. Due to the difficulty of applying the spin-boson model to multiple qubit 
systems, we take another approach. Instead of treating the bath explicitly, we follow the stochastic 
Liouville approach of the HSR model, and consider an effective Hamiltonian that treats the effect 



21 



22, 231. To describe the 



of the bath as a set of classical fluctuating fields acting on the system 
dynamics of an array of qubits under the influence of an external control field and environmental 
noise, we consider a system of n qubits, and start from a Hamiltonian with time independent and 
time dependent parts. The general Hamiltonian of the qubit system can be written as 
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H(t) = H + h(t) 

2»-l (!) 

= • %(';•>;,•,. 

i,j=0 

where cj and Cj are the creation and annihilation operators for the i-th state of the 2 n basis set. 
The time independent part Ho describes the interactions between qubits, while the time dependent 
part h(t) describes the fluctuations of the interactions due to the coupling to the environment. For 
simplicity, we assume throughout this work that the external control fields are switched on and 
off instantaneously, and the interactions introduced by the external control fields are constant in 
time; this corresponds to a rectangular pulse. More realistic pulse shapes can be incorporated into 
our treatment without too much additional work. In addition, a sequence of different rectangular 
pulses can be divided into time periods with a constant external field in each of them, and then 
treated separately using a different time independent Ho for each time period. Considering only 
constant external control fields does not affect the generality of this model. 

The time dependent part of the Hamiltonian describes the influences of the environment via 
fluctuations of the system energy. This term may include fluctuations from many different origins, 
such as the fluctuations of imperfect control fields, the fluctuations induced by the bath on the 
qubit excitation energy, the off-diagonal matrix element, and the inter-qubit interactions. Following 



Haken and Strobl 



211 ] . we consider the fluctuations as random Gaussian Markov processes with 



zero mean and 5-function correlation times: 



<M*)> = (2) 

(M*)MO> = Rij;krS{t-t>) 

Here the brackets () represents the thermal average over all bath degrees of freedom, and the time 
independent correlation matrix element Rij-ki is a real number describing the correlations between 
hij(t) and hki(t'). All Rij-ki elements form a 2 2n -dimensional correlation matrix R. In addition, 
we have the following symmetry property of R: 



Rij;kl — Rji;kl — Rij;lk — Rji;lk — Rkl;ij (3) 

The value of Rij-ki depends on the strength of the coupling to the environment; therefore, it is a 
measure of the noisiness of the environment. The 5- function correlation in time corresponds to 
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the limit of fast bath modulations, which assumes that the relaxation time of the bath is much 
larger than the characteristic time of the system. Therefore, this model should be valid at high 
temperature limit. Also note that although the effect of the temperature can be included by 
considering temperature dependent correlation matrix elements, there is no explicit temperature 
dependence in this model. We will discuss the consequences of this assumption and the applicability 
of this model in more detail in Section IVl 

The time independent part of the Hamiltonian Ho and the correlation matrix R determine the 
dynamics of the system. The values of Ho and R depend on the setup of the physical system, the 
various types of noise considered, and the nature of the bath. Note that in the HSR model, the 
system is limited to the one exciton subspace, and the matrix Ho and R can be obtained directly. 
However, in our n-qubit system, all 2 n states must be considered, and how to determine Ho and 
R is less obvious. In the following sections, we provide explicit examples of Ho and R for systems 
describing quantum teleportation and generic quantum gates. Generalization of the procedure to 
determine Ho and R for a general n-qubit system should be straightforward. Throughout this 
section we will only use the generic forms of Ho and R to derive the equation of motion that 
describes the time evolution of the n-qubit system under the influence of noise. 

The dynamics of the system is described by the stochastic Liouville equation (H = 1) 

p(t) = -i\K(t),P(t)], 

where p(t) is the density matrix of the system at time t. Using the statistical properties of h(t) 
[Eq. ((21)] and the symmetry property of the correlation functions [Eq. we can compute the 

exact equation of motion for the averaged density matrix elements of the system by applying the 
second order generalized cumulant expansion method to average over all fluctuations. The result 
we obtain is in a simple form: 

ftPafl = -i^HajPjfl + i^PajHjp 

3 j (4) 

-\"*^Rlk;kpPal ~ \"^^Rlk;kaf>W + '^ J Rpl;kaPkU 
A . / A . / A ■ - / 

where all the summations are over all 2 n state indices. In addition, we have defined the averaged 
density matrix of the system, p(t) = {p(t)). In Eq. ©, the dynamics of the averaged density 
matrix can be separated into a coherent part, due to Ho, and a incoherent part, due to R. The 
dissipation of the system is governed by incoherent dynamics, and is related to the elements of 
the fluctuation correlation matrix R. Notice that the form of Eq. © is similar to the form of 
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the widely used Redfield equation, with the relaxation matrix elements given by the corresponding 
Rij-ki terms in the equation j^j]. 

Eq. Q forms a system of 2 2n linear ordinary differential equations (ODE). Given the values 
of Hij and Rij-ku the ODE system can be solved efficiently to yield the time dependent averaged 
density matrix pit). In fact, in most one qubit and two qubit systems, the equations can be solved 
analytically, and the analytical formula for p(t) can be obtained. In general, we can calculate Ho 
and R from the Hamiltonian of the system and the correlations between fluctuations introduced 
by the environment. Once we have Ho and R, it is then trivial to solve Eq. © to yield a p(t) 
that fully describes the dynamics of the n-qubit system. This procedure is straightforward, and 
can be used to study the effect of noise in complex quantum computations. We demonstrate the 
applications of this model on the study of the effect of noise on quantum teleportation and general 
quantum two qubit gates in the next two sections. 



III. DISSIPATION IN QUANTUM TELEPORTATION 

Since first proposed by Bennett et al. in 1993 j^J, the concept of "quantum teleportation" 
has received much attention. By exploiting the entangled nature of an Einstein-Podolsky-Rosen 
(EPR) pair, a sender can transmit the quantum state of a qubit to a receiver, without physically 
transferring the qubit through space. In this section, we will apply our stochastic Liouville approach 
to study the effect of noise on quantum teleportation. 



A. Quantum teleportation 

We first consider the ideal scenario of teleporting one qubit from Alice to Bob. Suppose Alice 
and Bob share a EPR pair, labeled as qubit a and b, emitted from an EPR pair source, and Alice 
wants to teleport qubit c in state = cq\0) + ci|l) to Bob. The EPR pair source emits two 
entangled qubits in one of the four Bell states at time t = 0, and then the two qubits are sent 
through separate channels C a and to Alice and Bob, respectively. After receiving qubit a, 
Alice performs a Bell-state measurement on her qubits (a and c), and sends the outcome of her 
measurement to Bob through a classical channel. Alice's measurement projects qubit b onto one 
of the four corresponding states, i.e. I • (c |0) b + ci|l)&), a z ■ (c |0) b + ci|l) 6 ), a x ■ (c |0) b + ci|l)&), 
and io y ■ (co|0}& + ci | . Bob then applies the corresponding inverse transformation (I, a z , a x , 
and — ia y , respectively) to recover his qubit in the state 
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In practice, errors can happen during the quantum teleportation from several origins: (1) the 
noise in the quantum channels C a and C&, (2) the degradation of qubit c after the preparation, 
(3) the imperfect Bell-state measurement performed by Alice, (4) the further degradation of qubit 
b when transmitting the result of Bell-state measurement through the classical channel, (5) the 
imperfect unitary transformations performed by Bob. Here, we only consider the first situation 
where channel C a and Cb are noisy, and focus on the errors due to the degradation of entanglement. 
We assume all other operations are done perfectly. This means that result obtained in the following 
represents a lower bound on the errors in the quantum teleportation. 



B. Effect of noise on a pair of entangled qubits 

To study the degradation of a pair of entangled qubits, we consider the effective Hamiltonian 
describing two uncorrelated qubits a and b: 



H — H a + Hb 

= £ e n (t)-at ) + Y, J n(t)-4 U) ( ) 

n=a,b n=a,b ^ ' 

= i £ n + &n(t)] • 4 U) + E t J « + 5J "®\ ■ ^ 

n=a,b n=a,b 

where and o4 , n = a,b are Pauli spin operators on qubit a and b; 2e a (2e&) is the averaged 
energy splitting between the |0) and |1) states of qubit a (6); J a (Jf,) is the averaged off-diagonal 
matrix element for qubit a (6); 5e a (t) (Seb(t)) is the time-dependent fluctuating part of the diagonal 
energy for qubit a (6); 5J a (t) (5Jb(t)) is the time-dependent fluctuating part of the off-diagonal 
matrix element for qubit a (b). Following the assumption made in Section [HJ we regard Se n (t) and 
SJ n (t), n = o, b as Gaussian Markov processes fully described by their first two moments: 



(6e n (t)) = {SJ n (t)) = 0, 
(5e n (t)5e m (t')) = ^ • S nm S(t - t'), 

(6) 

(5J n (t)5J m (t')} = ~/?-6 nm 6(t-t'), 
(5e n (t)5J m (t')) = 0, 

where 7$ (7$ ) describes the strength of the diagonal energy fluctuations of qubit a (b); 7" (7J ) 
describes the strength of the off-diagonal matrix element fluctuations of qubit a (b). Clearly, 7$ 
and 7q are related to the system-bath interactions involving a z system operators, and 7° and 7J 
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are related to the interactions involv ing a x system operators. These phenomenological parameters 
can be estimated experimentally j^.!^. Notice that we treat the correlation between qubit a and 
b independently, because in quantum teleportation, the two EPR qubits are sent through different 
channels to two distantly separated places, thus the two qubits are coupled to distinct baths. In 
addition, we assume the diagonal and off-diagonal fluctuations are not correlated. 

To simplify our computations, we choose to study the dynamics of the system in the Bell-state 
basis. The four Bell states are defined as 

^(|0>«|0> 6 + |1)«|1) 6 ), 



\Bi) 

\B 2 ) = ^(|0)„|0) 6 -|l) a |l) 6 ), 

\Bz) = ^(|0) a |l) 6 + |l) a |0} 6 ), 

\BA = -4no)„ii) 



^VI"/a|-/6 - |l)a|0)ft), 

where subscript a, b labels the state of different qubits. For convenience, hereafter we will use the 
notation that use the first digit to denote the state of qubit a, and the second digit to denote the 
state of qubit b, i.e. |l) a |l){, = |11). The Hamiltonian for the two qubit system [Eq. (J5J)] in the 
Bell-state basis is 



H 



e a + e b + h 12 (t) J a + J b + h 13 (t) 

e a + e b + h 2 i{t) J b - J a + h&(t) 

J a + J b + h 31 (t) e a -e b + h 34 (t) 

J b - J a + h 42 (t) e a -£ b + h 43 (t) 



(7) 



where the nonzero transformed time-dependent matrix elements are: 



hi 2 (t) = h 21 {t) = 5e a {t) + 5e b (t), 
hi 3 (t) = h 3l (t) = 5J a {t) + 5J b (t), 

M*) = M*) = sj b (t)-sj a (t), 

h u (t) = h A3 (t) = 5e a (t) - 5e b (t). 

From Eq. © and Eq. (JHJ), we can easily compute the correlation matrix R of the system. In 
this case, R has only 32 nonzero elements that can be represented by the following 6 irreducible 
elements: 
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-Rl2;12 


= 7o a + 7o' 


-Rl2;34 


= 7o a -7o & 


-Rl3;13 


= 7? + 7i 


-^13;23 


= 7? " 7? 


-??24;24 


= 7? + ll 


-^34;34 


= 7o a + 7o 6 



(9) 



Other nonzero elements of R can be obtained using the symmetry property of R (Eq. ©)• 
Plugging the correlation matrix elements Eq. ([?))] and the time-independent Hamiltonian matrix 
elements [Eq. Q] into Eq. @, we obtain the equation of motion for the averaged density matrix 
of the system, p(t). 

In the limit of zero averaged Hamiltonian matrix elements, e n = J n = 0, the equation of motion 
for the diagonal density matrix elements are decoupled from those for the off-diagonal density 
matrix elements. Therefore, the dynamics of a system initially in one of the four Bell states 
(i.e. the initial density matrix has only non-zero diagonal elements) can be fully described by the 
equations for the diagonal density matrix elements: 



£ t p n (t) = T ■ [p 22 (t) - p n (t)\ + Tx • [p 33 (t) - pu(t)} , 

i t P22(t) = To ■ [pil(t) - P2l{t)\ + T l ■ [p U (t) - P22(t)} , 

&Paa{t) = T ■ [p4A(t) - paait)] + T x ■ [p n (t) - p 33 (t)} , 

f t pu{t) = To • [p33(t) - pu(t)\ + r x • [p 22 {t) - p U {t)} , 

where we have defined Tq = (7$ + 7o)i an( i = (7" + These equations have the form of 
a system of kinetic equations involving four states, and, clearly, To and I?i have the meaning of 
the degradation rate constants. The symmetric form of Eq. (|1U(I suggests that all four states are 
equivalent dynamically, hence we expect the degradation rates of the systems initially in any of the 
four Bell states are equal. In this limit, the results of the teleportation based on different Bell-state 
channels are the same. Later we will show that this is only true when e n = J n = and the two 
qubits are coupled to distinct baths. 

Eq. (|1U|) also shows that a system of two qubits initially in one of the maximumly entangled 
states degrades into a statistical mixture of the four Bell states. Assuming that the system is 
initially in the state \Bi) and stays in the noisy quantum channels for a time period t, the density 
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matrix for the entangled qubits Alice and Bob obtained can be represented as the statistical mixture 
p(t) = hl(t) ■ + m{t) • \B 2 )(B 2 \ + p 33 (t) ■ \B 3 )(B 3 \ + p M (t) ■ |B 4 )(5 4 |. (11) 

The populations can be obtained by solving Eq. (|T0"|) with the initial condition po = \Bi)(B\\: 



Pu(t) = 


1 
1 


+ 


¥ 


-2T t 


+ 


¥ 


-2rit 


+ 




-2(r +ri)t 


P22(t) = 


1 

4 




¥ 


-2Y t 


+ 


¥ 


-2Tit 






-2(r +ri)t 


P33(t) = 


1 
1 


+ 


¥ 


-2T a t 




¥ 


-2Txt 






-2(r +r!)i 


P4i(t) = 


1 
1 




¥ 


-2Y t 




¥ 


-2r x t 


+ 


¥ 


-2(r +ri)t 



(12) 



From Eq. (|12|) . the fidelity of the entangled pair, defined as the overlap between the initial 
density matrix po and the density matrix at time t, can be calculated: 

F e {t) = Trpop(t) = \ + \e-^ + \e~^ + \e~^^ . (13) 

Eq. (fT3)) shows that when F and Ti are both non-zero, the fidelity -F e (oo) = \ in the long time 
limit. When either To or T\ is zero, -F e (oo) = ^. This result indicates that if we can somehow 
transform the system and minimize either the diagonal energy fluctuations or the off-diagonal 
matrix element fluctuations, the original quantum state can be better preserved. In addition, Eq. 
1)1 3|) can be used to compute a critical time scale beyond which the degraded entanglement can 
not be purified by any entanglement purification method . The fidelity required by a successful 
entanglement purification process, F e {t) > 0.5, corresponds to a critical time t c where F e (t c ) = ^. 
For any high-fidelity quantum teleportation to be possible, the EPR pair should not be allowed 
to stay in the noisy channels for a time period longer than t c . t c also defines the critical distance 
for possible high-fidelity quantum teleportation, given the noise of the channels described by the 
parameters Tq and T\. 



C. Outcome of teleportation 

Now we can use the result in the previous section to study the outcome of teleporting a qubit 
c in state \ip) = co|0) + ci|l) from Alice to Bob. We assume the traveling time that the EPR pair 
spends in the noisy channels is t, and the averaged energy e n and off-diagonal matrix elements J n 
for both qubits are very small so that the limit of e n = J n = 0, n = a, b can be applied. After 
receiving the degraded EPR pair described by Eq. ()11)) . Alice and Bob then perform the Bell-state 
measurement and corresponding unitary transformation to complete the teleportation. Assuming 
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that all measurements and unitary transformations are carried out perfectly and do not introduce 
more error, the teleportation outcome that Bob obtains is 



At) 



1 + 1 (| CQ |2 _ | Cl |2) . e -2r lf £0£l+£5a . e -2T t + £o£^£5£l . ^(Fo+^t 

c cl+c* Cl _ e -2T t + cgci-cpc* _ e _ 2 (r +r 1 )t 1 + I(| Cl |2 _ | Cq |2) . g-2I\t 



(14) 

This result is similar to the result for the dissipation of a two level system in the HSR model 
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. Notice that the decoherence depends on the total diagonal fluctuations, Tq = 7q +7oi an d 
the population relaxation depends on the total off-diagonal fluctuations, Ti = 7J+7J. Clearly, noise 
in both channels affect the teleportation outcome additively. In fact, the outcome is exactly the 
same as if the teleported qubit is transfered physically from Alice to Bob through the noisy channel 
C a and Cj), although the qubit Bob receives has never traveled through channel C a physically. 
The fidelity of teleportation as a function of the traveling time t is 

F tele (t) = \ + ifcgd + coc^e" 2 ^ + i(|c | 2 - | Cl | 2 ) 2 e- 2r ^ - \^ Cl - c^V 2 ^*. (15) 

The fidelity of teleportation decreases monotonically from 1 to | as the traveling time t increases. At 
the long time limit, the fidelity approaches 5, which means the result of the quantum teleportation 
is a half-half mixture of |0) and |1) states, i.e. information about is totally lost. This result is 
in agreement with recent studies on the effect of noise on quantum teleportation j^ . 

Eq. (|15|) provides a simple interpretation for the phenomenological parameter Fq and I\: Tq is 
the total decay rate for the real part of the coherence, and Ti is the total population relaxation 
rate. Recall that 7q , n = a, b is defined using the second moment of the diagonal energy fluctuation 
Se n (t), n = a,b (coupling involving a z ) and 7" , n = a, b is defined using the second moment of the 
off-diagonal matrix element fluctuation 8J n (t), n = a,b (coupling involving a x ). We see clearly 
the effects of different types of noise: the diagonal fluctuations introduce phase shifts that only 
affect the coherence of the qubit; the off-diagonal fluctuations introduce coupling between the two 
states and result in population transfer. Note that the decay of the imaginary part of the coherence 
depends on both diagonal and off-diagonal fluctuations. In the terminology of quantum computing, 
phase-shift errors are caused by the diagonal energy fluctuations, bit-flip errors are caused by the 
off-diagonal matrix element fluctuations, and the change in the a y component are due to both types 
of fluctuations. Previous studies on the dissipation of qubits using spin-boson types of Hamiltonian 



give similar results for the effects of different types of system-bath interactions 



10, 
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11,03. Our 
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model gives direct relationships between the phenomenological parameters describing the strength 
of the fluctuations and the dissipation rates. In addition, our model can take into account the 
effects of both types of fluctuations simultaneously, which is different from most error models used 
previously. 

D. Nonzero averaged matrix elements 

When the time-independent part of the Hamiltonian contains nonzero matrix elements, i.e. 
e n ^ or J n 7^ 0, n = a, b, the exact analytical expression for p(t) is not generally available. 
In addition, the effect of diagonal energy fluctuations no longer can be clearly distinguished from 
the effect of off-diagonal matrix element fluctuations, both population relaxation and decoherence 
depend on 7q and 7", n = a,b . More importantly, the four Bell states no longer decay at 
the same rate, and we can see the effect of the coherent dynamics depending on the value of 
the averaged energy and off-diagonal matrix elements. In the weakly-damped regime where the 
averaged Hamiltonian matrix elements are larger than the strength of the noise, the dynamics of a 
pair of entangled qubits exhibits coherent oscillating behavior. These oscillations can lead to errors 
of the quantum teleportation. Figure ^ shows the fidelity of the four Bell states as a function of 
traveling time at e a = = 1, J a = Jb = 0.5, 7o = 7o = 0.1, and 7i = 7i = 0.1. The different 
oscillating behavior of the Bell states can be understood by considering the time-independent part 
of the Hamiltonian. From Eq. ©, all the nonzero time-independent matrix elements are 

(-B1IH0I-B2) = (-B2IH0I-B1) = E a + Eb, 

(-Bi|H |-B 3 ) = (B 3 \H \Bi) = J a + Jb, 

(B2IH0I-B4) = (-B4IH0I-B2) = Jb-Ja, 

(-B3IH0I-B4) = (.B4IH0I-B3) = e a — Eb- 

These matrix elements govern the coherent transition between the Bell states, and result in the 
oscillating behavior of the dynamics. In Fig. ^ the fidelity of the I-B4) state decays monotonically 
as t increases, because both matrix elements couple this state to the other states, Jf, — J a and 
£a — £b> are zero for the parameters used. This also explains why the fidelity of the I-B4} state 
provides an upper bound on the fidelity of other Bell states in Fig. ^ The state that is coupled 
most weakly to other states decay most slowly. 

In the regime where the averaged Hamiltonian matrix elements are smaller than the strength 
of the noise, the system is overdamped and no oscillating behavior can be observed. Figure 
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shows the fidelity of the four Bell states at e a = Sb = 0.1, J a = Jb = 0.05, 7o = 7o = 0.1, and 
7? = 72 = 0.1. In this regime, all Bell states degrade monotonically as the traveling time increases. 

The fidelity of the EPR pair used in the quantum teleportation is directly related to the 
fidelity of teleportation. Therefore, the above discussion can be directly applied to the fidelity of 
teleportation performed using different Bell states. When e n ^ or J n ^ 0, n = a, b , the fidelity 
of the teleportation behaves differently when different Bell states are used. To achieve the best 
result for the teleportation, we have to choose the Bell state that is coupled most weakly to other 
states. In general, e n > 0, n = a, b and J a and Jb have the same sign, thus I-B4) state will have the 
weakest coupling. The singlet \B$) state is the preferred EPR state for the quantum teleportation. 

[Fig. [1] 

[Fig. m 

E. Effect of collective bath 

We have studied the dissipation of two entangled qubits each coupled to a distinct bath, which 
is the typical situation relevant for the quantum teleportation. Another interesting case is when 
the two qubits are coupled to a common bath. In this case, we use the Hamiltonian of Eq. Q; the 
difference in the state of the bath is reflected by different correlation functions for the stochastic 
processes. When the two qubits are coupled to a common bath, the first two moments can be 
represented as 

(5e n (t)} = (5J n (t)) = 0, 
(5e n (t)5e m {t')) = lQ .5{t-t'), 

(16) 

(5J n (t)5J m (t')} = jx-S(t-lf), 
(5e n (t)5J m (t')) = 0, 

where 70 describes the strength of the diagonal fluctuations; 71 describes the strength of the off- 
diagonal fluctuations. Note that because the qubits are coupled to a common bath, the fluctuations 
on different qubits are correlated. From Eq. (JHJ) and Eq. (|16|). we can derive the correlation matrix 
R for the system in the Bell-state basis. In this collective bath limit, R has only 8 nonzero elements 
that can be represented by the following 2 irreducible elements: 

-Rl2;12 = 47o, 
#13;13 = 47i. 
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Using Eq. (|17[). we can derive the equation of motion for the dynamics of two qubits coupled 
to a common bath. In the limit of zero averaged Hamiltonian matrix elements [e n = J n = 0, 
n = a, 6), we obtain a simple result for the populations in the four Bell states: 

!/5 u (t) = 4 7o • " + 4 7i • feW - Pu(*)] , 

f t ~p 22 (t) = 4 7 o-[pii(t)-p 22 a)], „ . 

|p 33 (*) = 4 7l -[ / 5 11 (t)-p 3 3(t)], 

&m(t) = o. 

Eq. (|1H|) describes the dynamics for a system of two qubits coupled to a common bath in the 
Bell-state basis. Interestingly, the population in the \B±) state, p^it), is invariant in time. In 
addition, when only diagonal energy fluctuations exist (71 = 0), the population in the I-B3) state 
is also invariant; when only off-diagonal matrix element fluctuations exist (71 = 0), the population 
in the I-B2) state is invariant. Compared to the result of two qubits coupled to distinct baths [see 
Ea. (|lU|) ]. Eq. (|18|) shows that the fluctuations interfere constructively for the \Bi) state leading to 
a faster decay rate, but destructively for the I-B4) state. This result can be understood easily in 
our stochastic model. In our model, the effect of environment on the system is represented by a 
fluctuating field, and the interaction Hamiltonian for the two qubits is Hint = ■V a (t)+af' -Vb(t) 
(i = x, z; a and b are labels for different qubits). When the two qubits are coupled to a common 
bath, V a (t) = Vb(t), we can factorize the interaction into the form Hint = + ) • V a (t). 

Therefore, any state \4>) that satisfies (<^|c7- a ^ + <r- \4>) = does not interact with the fluctuating 
field, and is invariant to the noise. We can see that (B^a^ -\-a^\B^) = and (B^a^ +cr^ I-B4) = 
0, thus both the I-B3) and I-B4) states are not affected by phase-shifting noise; (B2\cr^ \B2) = 
and (B^\ai a ^ + &x \B^) = 0, thus both the \B2) and ^4) states are not affected by bit-flipping 



noise. This effect of the collective bath has been verified experimentally [29.1, and studied in 
theoretical works related to the ideas of "quantum error-avoiding codes" 30, 31J and "decoherence- 

nn 

free subspaces" 32, 33]. Duan et al. have shown similar result using a Hamiltonian that explicitly 
includes the linear coupling terms between the system and the boson bath jsil . I34I . The agreement 
indicates that our simple stochastic model can handle both the independent and the collective bath 
properly. 

Recently, Kumar and Pandey have studied the effect of noise on quantum teleportation [2^ . 
They applied two different models, a stochastic model and a spin-boson type model, to this problem, 
and studied the relative teleportation efficiencies of the Bell states. Their main result is that for 
the simple stochastic model, the four Bell states are equivalent, but for the second model in which 
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the effect of environment is considered explicitly, the \B±) state is least affected by the noise. We 
obtain a similar conclusion using the stochastic Liouville equation approach. Based on our result, 
we understand that the \B&) state is the least affected state because of the assumption of a collective 
bath, not because the effect of bath is considered microscopically. Like spin-boson type models, 
a simple stochastic model when treated correctly can provide the same result, and gives a simple 
picture for the effect of a collective bath versus a localized bath. 



IV. ERRORS IN QUANTUM CNOT GATE 

Qubits and quantum gates are the basic elements of quantum computing. A quantum circuit 
that performs a particular quantum operation can be expressed as a composition of elementary 
quantum gates [351 ] . In fact, quantum circuits can be constructed using one- and two-qubit gates as 
basic building blocks. For example, the quantum CNOT gate together with all one-qubit quantum 
gates form such a set of universal quantum gates _.36<|. In reality, quantum computations are 
performed by subjecting an array of qubits under a sequence of control fields that control the 
Hamiltonian of the qubit system and result in specific quantum gate operations. Therefore, we 
consider the process of quantum computation as preparing the qubit system in the initial state, 
then performing programmed control fields on the qubits in a sequence of time steps, and finally 
measuring the output in the working basis. 

To understand the effect of noise on general quantum computations and help the implementa- 
tion of quantum computers, we need a model that can be used to describe the decoherence and 
population relaxation for a system of qubits subjected to external control fields. The decoherence 
and gate perf ormance of a CNOT gate on various types of physical realizations have been studied 
in Refs. 14, Irsl . ItgI . IitI . In particular, Thorwart and Hnggi investigated the decoherence 

and dissipation for a generic CNOT gate operation using the numerical ab initio technique of the 
quasiadiabatic-propagator path integral (QUAPI). They demonstrated that this numerical method 
is capable of describing the full time-resolved dynamics of the two-qubit system in the presence of 
noise. To our knowledge, so far, the QUAPI method is the most sophisticated method that has 
been applied to study the decoherence during a CNOT gate operation. In this section, we apply 
the stochastic Liouville equation approach to study the same generic CNOT operation investigated 
by Thorwart and Hnggi, and show that our model yields similar results. In general, our model is 
easier to extend to many qubit systems than the QUAPI method, and can incorporate the effects 
of noise from different sources at the same time. 
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A. A generic model for 2-qubit quantum gates 

In a physical system, a quantum gate can be expressed by a Hamiltonian with terms representing 
the control fields that result in the gate operation. Consider a elementary step in a quantum gate 
operation where the control Hamiltonian is switched on, a generic Hamiltonian describing the 
constant external fields and the time-dependent fluctuations (noise) for a two-qubit system can be 
written as 



H(t) = Yl [ £ « + «M*)] • 4 n) + E [ J « + 6J ^ ■ a * ] 

n=a,b n=a,b 

= H + h(t) 

where the two qubits are labeled as qubit a and qubit b; the first two terms comprise the Hamilto- 
nian for two non- interacting qubits considered in Eq. the last term represents the inter-qubit 
interaction with o±^ = (a^ =F toy )> n = a,b; g and ^(i) are the time-independent and time- 
dependent fluctuating part of the inter-qubit coupling. The controllable fields are represented by 
the values of s n , J n , n = a,b, and g. Quantum gates can be implemented by switching these fields 
on and off in a controlled manner. Notice that the XY type of coupling is adopted in our model 
Hamiltonian. This interaction is just an illustrative example, and does not account for all the 
possible interactions in a specific realization of solid-state devices. The real form of the inter-qubit 
interaction term depends on the controllable interactions available for each individual physical im- 
plementation. Nevertheless, our model can handle the other types of interactions as well, and we 
expect that the model Hamiltonian we use here can reproduce the same general physical behavior 
as other two-qubit Hamiltonians. 

From Eq. (|19|) . we can write down the time-independent part of the Hamiltonian in the standard 
basis {|00), |01), |10),|11)}: 



Hn 



£a + £6 


Jb 


Ja 





Jb 




9 


Ja 


Ja 


9 


£fc — £a 


Jb 





Ja 


Jb 


—e a — £b 



(20) 



and the time-dependent part of the Hamiltonian is 
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h(t) 



<5e a (i) + fe 6 (t) 
SJ b (t) 

5Ja(t) 



SJ b (t) 
5e a (t) - 5e b (t) 



Sg(t) 
8J a (t) 



Se b (t) - 5e a (t) 
SJ b (t) 



SJ a (t) 
Sg(t) 



SJa(t) 
5J b (t) 

5e a (t) - 8e b (t) 







(21) 







Furthermore, we assume the two qubits are close to each other in space, therefore, we consider 
the correlation functions suitable for two qubits coupled to a common bath. Again, we assume the 
fluctuations have zero mean and 5- function correlation times. The nonzero second moments are 



where 70 describes the strength of the diagonal energy fluctuations; 71 describes the strength of the 
off-diagonal matrix element fluctuations; 72 describes the strength of the fluctuations of the inter- 
qubit interactions. As we have shown in the previous section, these phenomenological parameters 
are related to the kinetic rate of each separate dissipative process, and can be easily measured 
experimentally. Also note that we directly include the inter-qubit coupling fluctuations, which 
corresponds to two-qubit flip-flop errors that are difficult to treat in the microscopic spin-boson 
type Hamiltonians. 

Eq. ()22|) can be used to compute the elements of the correlation matrix R. Using R together 
with the averaged Hamiltonian matrix elements in Eq. (|2Uj) , we can obtain the equation of motion 
describing the dynamics of the two-qubit system subjected to arbitrary one- and two-qubit control 
fields. As a result, we can study the dissipative dynamics of the qubit system during arbitrary gate 
operations. Although we only consider an operation done by a set of constant external fields, the 
behavior of more complicated gates that involve more than one step can be studied by combining 
the result for each elementary operations. In our model, the results for a set of universal quantum 
gates can be assembled to compute the results for a general quantum circuit. 



(5e n (t)8e m (t')) = l0 -5(t-t') 
(5J n (t)5J m (t')) = Tl • 5(t - i') 
(5g(t)5g(t>)) = l2 -5{t-t>) 



(22) 



B. The quantum CNOT gate 



The quantum CNOT gate plays a central role in the quantum computation, because, as we 
noted above, the set of all one-qubit gates together with the CNOT gate is universal j^. In the 
standard basis {|00), |01), |10),|11)}, the ideal CNOT gate is represented as 
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10 
10 

jTideal 
u CNOT 

1 
10 

This gate operates on two qubits, and inverts the state of the second qubit if the first qubit is in 
the state |1). The CNOT gate cannot be constructed in one step using our model Hamiltonian. 
Instead, we must construct the CNOT gate using multiple elementary one- and two-qubit gates. 
To begin with, we define the following one-qubit rotations on qubit a and b: 



Unz(a) 



;«>) 

ea z , n = a,t 



» 



(a 



n 



a,b, 



and the two-qubit operation: 



. , (a) (b) (a) (6), 

Uj(a) = e ta ^+ a - +a - a + \ 



All these operations can be easily implemented using our model Hamiltonian [Eq. (|19|)] (with all 
control fields set to zero initially): U nz (a), n = a,b, can be done by switching on e n = —e$ ■ Sign(a) 
for a time period of r = U nx (a), n = a,b, can be done by switching on J n = —Jo • Sign(a) 
for a time period of r = Uj(a) can be done by switching on g = —go ■ Sign(a) for a time 
period of r = where the sign function Sign(a) returns —1 when a < 0, and 1 when a > 0. 
Using the corresponding averaged Hamiltonian Ho for each operations and the correlation matrix 
presented in the previous section, the equation of motion describing the dynamics of the two-qubit 
system subjected to any of these operations can be easily obtained. Actually, for arbitrary initial 
conditions, the analytical solution for the time-dependent two-qubit density matrix pit) during 
U nx {ot), U nz (a), n = a,b, and Uj(a) operations are available in the Laplace domain, and can be 
used to study arbitrary quantum circuits composed by these three elementary operations. 

The CNOT gate can be expressed by the following sequence of one- and two-qubit gate opera- 
tions j| : 



UcNOT = UUlPU^PU-^Uji^Ua^^U^UU^PaA^)- (23) 
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Table (QJ) lists the required control fields and time span to implement each step using our model 
Hamiltonian. In Table (QJ, we use £o, ^0) an d go to denote the strength of the controllable 
single-qubit bias, intra-qubit coupling, and inter-qubit XY interaction, respectively. In addition, 
we assume that the controllable field strengths and noise (defined by parameters 70, 71, and 
72 as mentioned in the previous section) for the two qubits are identical. The value of these 
parameters should depend on the specific physical realization of the qubit systems. The total 
time required to perform the CNOT gate is T cnot = tt/2eo + vr/Jo + it /go. For a typical energy 
scale of 1 meV (suitable for quantum dot qubits), the operation time is on the picosecond time scale. 



[Table |U 

Using the parameters listed in Table (QJ), we can calculate the time-dependent two-qubit density 
matrix pit) during CNOT operations under different noise conditions defined by 70, 71, and 72. 
Figure 01 shows the time-resolved CNOT operation for two qubits initially in the |11) state. We set 
the strengths of the control fields equal to 1, i.e. Eo = Jo = go = 1. The ideal operation (solid line) 
starts at population 1 in the |11) state, and ends its total population in the 1 10) state, showing a 
successful CNOT operation. Three different noisy operations are shown in Fig. |HJ (1) operation 
with the strength of the diagonal energy fluctuations 70 = 0.05 (dashed line), (2) operation with the 
strength of the off-diagonal matrix element fluctuations 71 = 0.05 (dash-dotted line), (3) operation 
with the strength of the inter-qubit coupling fluctuations 72 = 0.05 (dotted line). The effect of 
noise on the CNOT operation can be clearly seen. In previous work, Thorwart and Hnggi derived 
the same time-resolved CNOT operation result Our result is very close to their numerical ab 
initio QUAPI result. The agreement between our time-resolved result to the QUAPI result gives 
us confidence that our model captures the correct physics. 

We use the gate fidelity and gate purity to characterize the performance of the CNOT gate. 
Other gate quantifiers including the quantum degree and entanglement capability are also calcu- 
lated |38j], but we do not show the results here because they follow the same trend as the gate 
fidelity and gate purity. In our formalism, the density matrix for the two qubits after the noisy 
CNOT operation, p(T cno t)=UcNOTPoUQ NOT , can be calculated for any initial density matrix po- 
Following Thorwart and Hnggi , we average the gate fidelity and gate purity over 16 initial states 
to account for the general performance of the CNOT gate. The 16 unentangled input states \i/Jq), 
i,j = 1,2,3,4 are defined as = |&}«® \4>j)b with \4>x) = |0>, |0 2 ) = |1>, \<j> 3 ) = (|0> + \l))/y/2, 
\4>s) = (|0) + \l))/y/2, 1 04) = (|0) + i\l))/y/2, and a, b denoting the state for different qubits. These 
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states span the Hilbert space for the two-qubit operations, and should give a reasonable result for 
the averaged effect Q| . 

The gate fidelity is defined as the overlap between the ideal output and the output of the real 
gate operation. Using the 16 initial states, the averaged fidelity can be written as 

1 ; 

F =JqY1 (^outlPcNoMit)' 

where we have defined the ideal CNOT output \tp l Jut) = ^cnotI^o ) ^ anc ^ the output of the real 
CNOT operation P 1 c NO t = Ucnot\4 ) q){^q\Uq N q T - The gate fidelity is a measure of how close the 
real operation is compared to the ideal operation. For a perfect gate operation, the gate fidelity 
should be 1. 

Similarly, the averaged gate purity is defined as 

P =^E Tr ((^OT) 2 )- 

ij=l 

The gate purity quantifies the effect of decoherence. For a perfect gate operation, the gate purity 
should be 1. 



[Figure |3] 

C. Dependence on the noise strength 

The results of the averaged gate fidelity and gate purity as a function of the strength of each 
individual type of noise are shown in Fig. ^ For our generic study, we again set the strengths of 
all the control fields to 1, i.e. £0 = ^0=50 = 1- Clearly, different types of noise cause different 
amount of errors. However, they all follow the same trend. The deviations of the gate fidelity and 
gate purity from the ideal values, i.e. 1 — F and 1 — P, are sensitive to the strength of the noise, and 
saturate to 0.75 in the strong noise limit; the value 0.75 corresponds to a fully mixed state. In the 
weak noise regime, both 1 — F and 1 — P depend linearly on the noise strength, as expected lflllfil|. 
The proportionality constant in this case is ~ 10. In fact, the proportionality constant depends 
on the strengths of the control fields, and reflects the total operation time required to complete 
the CNOT gate operation. As the strength of the control field increases, the total operation time 
decreases, and the qubits have less time to undergo the dissipative processes, resulting in less 
degradation. To minimize the effect of noise, we need to reduce the proportionality constant, 
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therefore, we will want to operate the device at the highest control fields possible. However, the 
situation will be different if increasing the strengths of the control fields will also introduce more 
noise. We will explicitly discuss the effect of the control-field strength in the next subsection. 

From our results for £q = Jo = 9o = 1, to achieve the threshold accuracy of the 0.999 99 level 
needed for arbitrary long quantum computations [jjj], I40I 41 ] , one needs to keep the noise strength 



below the 10 -6 level. Assuming a characteristic energy scale of 1 meV, this value corresponds to 
a decoherence time 7" 1 in the /jls scale, which provides a serious challenge for experimentalists 
working on the realization of solid-state quantum computers. 



[Fig. m 

The linear dependence oil — F and 1 — P on the noise strengths also indicates that the effect 
of the same type of noise is additive in the weak noise regime. To study the additivity of different 
types of noise, we calculate the averaged CNOT gate fidelity when different types of noise coexist 
at the same time. We define the total error of the CNOT gate operation E as the deviation of the 
gate fidelity from the ideal value: 



#(7o,7i,72) = 1-^(70,71,72), (24) 

where we have explicitly expressed the total error £ as a function of the three different types of 
noise strengths:7o, 71, and 72. In Fig. 03 we show the errors of the CNOT gate operation where 
the different types of noise coexist, and compare them to the total errors obtained by adding up 
the errors caused by the individual type of noise. Clearly, for all four situations considered, these 
two lines collapse in the weak noise regime. The results indicate that errors caused by different 
types of noise are additive in the small noise regime. In other words, the following identity holds 
in the small noise regime: 

#(7o, 71, 72) = #(7o, 0, 0) + E{0, 7l , 0) + E{0, 0, 72 ). (25) 

Eq. 1(231) justifies previous studies where different types of system-bath interactions are treated 
independently lf|]. 



[Fig. E] 
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D. Dependence on the strength of the inter-qubit coupling 

In the solid-state implementations of qubit devices, the achievable inter-qubit coupling (go in 
our model) is usually weaker than other single-qubit interactions. Since the time required to finish 
a quantum gate operation is inverse proportional to the strength of the control field used, the weak 
interaction means long operation time, hence more errors. The two-qubit gate operations (Uj(a) 
in our model) is usually the bottleneck of quantum computation. In this section, we analyze the 
dependence of the quality of the quantum CNOT gate operation on the strength of the inter-qubit 
coupling go- 
lf the strength of the inter-qubit coupling go can be increased without introducing any extra 
disturbance on the system, then we expect operating the device in the strongest go achievable 
will give the best result. However, physically, applying a stronger field also means introducing 
stronger noise due to the imperfectness of the field. In our model, this means stronger fluctuations 
on the inter-qubit XY interaction. The extra noise can be expressed in the value of the 72 term. 
To incorporate this effect, we allow 72 to depend on the strength of the inter-qubit coupling 
go. Figure El shows the errors of the CNOT gate operation as a function of go at 70 = 0.001, 

71 = 0.001, Eo = 1, and Jo = 1. Three different noise strength dependences are shown: (1) 
constant 72 = 0.001 (solid curve), (2) linear 72 = 0.001 • (1 + go) (dashed curve), and (3) quadratic 

72 = 0.001 • (1 + <?o) (dash-dotted curve). The three curves show the same behavior in the 
small go regime, in which the operation takes too much time and the system is fully degraded. 
As the strength of the coupling go increases, the errors decrease due to the shorter operation 
time. When the strength of the coupling go approaches the strengths of other control fields 
(eo = Jo = 1 in this case), the three curves start to show different behavior. For both constant 
and linear 72, the errors generated by other operations (U nz (a) and U nx (a)) dominate the errors 
of the CNOT gate operation, therefore, increasing go gains nothing and the curve saturates. 
Our result for the constant 72 case is in agreement with the result obtained previously using 
the QUAPI method 15]. The situation is different when the strength of the noise depends on 
go quadratically. For this case, the errors start to increase after go > 1, because increasing 
the inter-qubit coupling go introduces stronger noise that cannot be compensated by shorter 
operation times. Therefore, in the quadratic case, there exists an optimal go for the gate operation. 



[Fig. © 
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V. LIMITATIONS AND POSSIBLE EXTENSIONS 

We have shown that the generalized HSR model is flexible for realistic physical devices. Appli- 
cations of this model on the effect of noise on the quantum teleportation and CNOT gate operation 
give us similar results compared to previous studies based on microscopic models. In this section, 
we will briefly discuss the limitations and possible extensions of this stochastic Liouville equation 
approach. 

A key step in the HSR model is to replace the microscopic system-bath interactions by stochastic 
processes. This procedure has permitted a full description of the dissipative dynamics of qubit 
systems and their response to the external fields. At the same time, we introduce phenomenological 
parameters to describe the strengths of fluctuations (70, 71 , and 72 in our model). These parameters 
have to be determined experimentally or computed using a separate microscopic model (22 , 42 , 4(J • 
Generally, 70, 71, and 72 should depend on temperature and increase as temperature increases. 
However, our model lacks explicit temperature dependence for these parameters, thus cannot be 
used to study the temperature dependence of the qubit dynamics. Fortunately, these parameters are 
directly related to physically measurable quantities, and can be easily determined by experiments. 
In our model, 70, 71, and 72 correspond to the decoherence rate, population relaxation rate, 
and inter-qubit flip-flopping rate, respectively; all of them can be measured by one- and two- 
qubit experiments. In addition, recent theoretical studies on the temperature dependence of the 
quality of quantum CNOT gate operation suggest that the temperature dependence of the gate 



performance is weak |l5l . 16], which is reasonable in the weak coupling regime and the temperature 
range relevant to solid-state qubit systems. 

The assumption of the fast modulation of the bath might be a more serious problem for the HSR 
model. The 5-function correlation time corresponds to an infinite fast decay of the bath correlations, 
which leads to incorrect short time dynamics and long time equilibrium populations. Palma et al. 
have studied the decoherence of a qubit and shown that the dynamics exhibits a "quiet" and a 
"quantum" regime at short times, and a "thermal" regime at long times The HSR model 

assumes that the bath relaxes infinitely fast, thus neglects the dynamics of the system before bath 
relaxation takes place. Although the HSR model cannot predict the short time dynamics correctly, 
we expect the physics for longer operations important for quantum computing are reasonably well 
captured. The delta function correlation can be replaced by an exponential function in time, and 
the extended model for a dichotomic process has been solved exactly without further assumptions 
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46j,|42t- 

It will be interesting to apply these extended models to quantum computations 
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and compare the results with the 5 correlation function results. 

The white noise assumption in the HSR model also corresponds to a bath with infinite tem- 
perature, therefore, the resulting equation of motion does not satisfy detailed balance at finite 
temperatures. As a consequence, the system always relaxes to equal populations regardless of the 
energy differences between the states. Extensions of the HSR model to solve this problem has been 
proposed in Ref. 4^ |. In quantum computing, we are mainly concerned about the dynamics of 
an unbiased qubit system, and even when a bias field is applied to the system to perform gate 
operations, the time period has to be short to avoid any population relaxation. Since we will only 
operate the quantum computer in the time scale that the population relaxation is negligible, we ex- 
pect the violation of the detailed balance condition will not cause serious problems for applications 
related to quantum computing. 

The stochastic representation for the dynamics of a quantum two-level system has been in- 



vestigated in Refs. 



481 ] and 49]. The correspondence between the phenomenological parameters 



describing the stochastic field (70 and 71 in this work) and the two-level system microscopic quan- 
tities are also studied. The stochastic approximation is found to be able to reproduce the results 
by Leggett et al. for the spin-boson model [lsl ]. Our results presented confirm this observation. 
In general, the stochastic Liouville equation approach presented in this work is applicable in the 
weak system-bath interaction limit relevant to quantum computations. 



VI. CONCLUSION 



In this work, we present a stochastic Liouville equation approach that provides a simple way 
to evaluate the effect of noise in quantum computations. This approach is generalized from the 
HSR model. Using an effective system Hamiltonian that includes the system-bath interactions 
as stochastic fluctuating terms with zero mean and delta function correlation times, we derived 
the exact equation of motion describing the dissipative dynamics for a system of n qubits. This 
generalized equation of motion is similar to the form of the widely used Redfield equation, with 
the relaxation matrix elements given by the corresponding correlation matrix elements. 

We then applied this model to study the dissipative dynamics of a system of two independent 
qubits that mimics the EPR pair used in the quantum teleportation. We showed that the phe- 
nomenological parameters used in our model, i.e. 70 and 71, correspond to the decoherence and 
population relaxation rate, respectively. To study the effect of noise on quantum teleportation, 
we calculated the fidelity of quantum teleportation. We found the effect of noise in the quantum 
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channels are additive, and the teleportation fidelity depends on the state of the teleported qubit. 
When the two EPR qubits are degenerate and have no intra-qubit coupling, the relative efficiencies 
of teleportation for the four Bell states are the same; otherwise, the singlet state I-B4) is the most 
efficient one. When the two qubits are coupled to the same bath (collective decoherence case), the 
\B\) state is superdecoherent, while the I-B4) state is decoherence-free. 

Furthermore, we studied a generic two-qubit Hamiltonian containing XY type inter-qubit inter- 
action. The dissipative dynamics of a set of one- and two-qubit quantum gates were studied, and 
the results were then combined to calculate the averaged gate fidelity and gate purity for the quan- 
tum CNOT gate operation. The dependence of the quality of the quantum CNOT gate operation 
on the noise strength and the strength of the inter-qubit coupling were investigated. We found that 
the quality of the CNOT gate operation is sensitive to the noise strength and the strengths of the 
control fields. In addition, the effect of noise is additive regardless of its origin. We compared our 
results to Thorwart and Hnggi's results obtained by the numerical ab initio QUAPI technique. In 
general, our results are in good agreement with those obtained by the numerical QUAPI method. 

We also discussed the limitations of the HSR type approach. The consequences due to the 
procedure of replacing the system-bath interactions by classical fluctuating fields and the assump- 
tion of the white noise were considered, and the possible extensions were noted. Generally, the 
application of HSR type model in the weak coupling regime that is relevant to quantum computing 
is justified. 

Finally, we emphasize that the model presented in this work can be used to study the dissipative 
dynamics of a many-qubit system with direct inter-qubit coupling, imperfectness of the control field, 
and other many-qubit effects. The power of this method is in its simple physical interpretation of 
the parameters and dissipative dynamics. In addition, because of the <5-function correlation time 
assumed in the model, there is no time-ordering problem for the propagator computed using Eq. 
P|. The resulting propagator satisfies complete positivity, therefore no additional time period has 
to be inserted between switching events, as will be necessary for methods based on the Bloch- 
Redfield formalism. As a result, propagators computed for simple one- and two-qubit gates can be 
directly assembled to study the dissipative dynamics of more complicated quantum circuits. Since 
this model can handle noise that affects multiple qubits at the same time, e.g. the fluctuations 
in the XY type inter-qubit interaction that results in flip-flop errors, it will be very interesting 
to apply this method to analyze the behavior of a quantum circuit implementing quantum error- 
correcting or fault-tolerant codes under the influence of various multiple-qubit noise. We also 
expect this method to be applied to evaluate the quality of quantum circuits under realistic device 
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conditions. Such theoretical studies will be useful for the design and implementation of quantum 
computers. 
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TABLE I: Parameters of the model Hamiltonians used to perform the CNOT gate in 7 steps. The required 
control fields and time span for each step are listed. Note that we only list the nonzero field parameters. 



No. 


Operation 


Control Fields 


Time 


1 


U bt (=?)U at {=£) 


e a = so, e b = e 




2 




g = -go 


^ = ^ + w. 


3 




Ja = Jq 




4 




g = go 


^ = ^ + Wo 


5 


Ub x (-n) 


Jb = Jo 


T 5 = T 4 + 2^ 


6 




£b = £o 


^ = ^ + ik 


7 


t4*(f) 


Jb = — Jo 
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Figure Captions: 

FIG. Fidelity as a function of the traveling time for the Bell states in the coherent 
regime: e a = e b = 1, J a = J b = 0.5, 7^ = 7^ = 0.1, and 7^ = 7^ = 0.1. 

FIG. EJ Fidelity as a function of the traveling time for the Bell states in the over-damped 
regime: e a = e b = 0.1, J a = J b = 0.05, 7o = To = 0.1, and 7? = 7l fe = 0.1. 

FIG. H3 Tim-resolved CNOT gate operation on the |11) input state. Shown are the popu- 
lations in the four basis states P«(£) = (ij\p{t)\ij) as a function of time. The strengths of all 
the fields are set to 1 in the calculation, i.e. £0 = Jo = 9o = 1, an d the corresponding time 
steps are defined in Table HJ). We show the results for four different CNOT gate operations: (1) 
ideal operation without any noise (solid line), (2) operation with the strength of the diagonal 
fluctuations 70 = 0.05 (dashed line), (3) operation with the strength of the off-diagonal fluctu- 
ations 71 = 0.05 (dash-dotted line), (4) operation with the strength of the inter-qubit coupling 
fluctuations 72 = 0.05 (dotted line). 

FIG. Dependence of the errors in the CNOT gate operation on the noise strength. The 
deviations of the gate fidelity (upper panel) and gate purity (lower panel) from the ideal values 
are shown, i.e. 1 — F and 1 — P. The effects of three types of noise are shown in both plots: (1) 
diagonal fluctuations represented by 70 (solid line), (2) off-diagonal fluctuations represented by 71 
(dashed line), (3) inter-qubit fluctuations represented by 72 (dash-dotted line). The control- field 
strengths are set to £q = Jq = go = 1. 
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FIG. 03 We show the error functions -E(7o, 71,72) of the CNOT gate operation in situations where 
the different types of noise coexist (solid lines). For each case, the corresponding total errors 
obtained by adding up the errors caused by the individual types of noise is also shown (dotted 
lines). Four different combinations are compared: upper-left: .E(r,r,0) vs. E(T,0,0)+E(0,T,0) 
(70 and 71); upper-right: E(0,T,T) vs. E(0, V, 0)+E(0, 0, T) (71 and 72); lower-left: E(T,0,T) vs. 
£(r,0,0)+£(0,0,r) (70 and 72); lower-right: E(T,T,T) vs. E(T, 0, 0)+E(0, T, 0)+E(0, 0, T) (all 
types of noise). The strengths of all the control fields are set to 1, i.e. £o = Jo = 9o = 1- We can 
clearly see that errors caused by different types of noise are additive in the small noise regime. 

FIG. Dependence of the errors in the CNOT gate operation on the strength of the 

inter-qubit coupling qq. Shown are the deviations of the gate fidelity from the ideal value for 
three types of 72: (i) constant 72 = 0.001 (solid curve), (ii) linear 72 = 0.001 • (1 + go) (dashed 
curve), and (iii) quadratic 72 = 0.001 • (1 + g{j) (dash-dotted curve). Other parameters are set to 
70 = 0.001, 71 = 0.001, e = 1, and J = 1. 
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